Enhanced Production of A and S(1385) Resonances 
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Yields of A(1230), E(1385) resonances produced in heavy ion collisions are studied within the 
framework of a kinetic master equation. The time evolution is driven by the process A ^^ Nn, 
E(1385) ^^ Att. We obtain resonance yield both below and above chemical equilibrium, depending 
on initial hadronization condition and separation of kinetic and chemical freeze-out. 

PACS numbers: 24.10.Pa, 25.75.-q, 25.75.Nq, 12.38.Mh 



Hadron resonances are produced copiously in the 
quark-gluon plasma (QGP) fireball break up into hadrons 
(hadronization, chemical freeze-out) e.g. at RHIC [l|, 0, 
BHi- Within the statistical hadronization model (SHM) 
approach [5|, Q , the initial yields are described by chem- 
ical fugacities T, and hadronization temperature T. The 
production of heavy resonances is suppressed exponen- 
tially in m/T. Once formed, resonances decay. If this 
occurs inside matter, detailed balance requires also pro- 
duction of resonances, called 'regeneration' and/or 'back- 
reaction'. 

If the chemical freeze-out occurs much earlier than 
thermal, the initially produced resonances are practi- 
cally invisible due to rescattering of decay products Q- 
The observed yield of resonances is fixed by the physical 
conditions prevailing at the final breakup of the fireball, 
at which time last scattering occurs, this is the 'kinetic 
freeze-out'. The present work addresses two questions: 

a) how observable resonance yield depends on the differ- 
ence between chemical freeze-out temperature (e.g. point 
of hadronization of QGP) and the kinetic freeze-out tem- 
perature; 

b) how this yield depends on the degree of initial chemi- 
cal non-equilibrium at hadronization. 

One can see this work as an effort to improve on the 
concept of chemical freeze-out for the case of resonances: 
given the relatively fast reactions their yield remains sen- 
sitive to the conditions prevailing between chemical and 
thermal freeze-out, even if this time is just Ifm/c. 

In order to describe evolution of the resonance abun- 
dance one can perform a microscopic transport simula- 
tion of the expanding system. In this approach the regen- 
eration of resonances was previously studied by Bleicher 
and collaborators 8, 9, 10]. There are many detailed fea- 
tures of particle interactions to resolve in a microscopic 
model description and thus it seems appropriate to sim- 
plify the situation. We study resonance decay and re- 
generation using the momentum integrated population 
master equations, and assuming hydrodynamic expan- 
sion inspired model of fireball dynamics with conserved 
entropy content. In all our considerations we presume 
that the yield of pions tt is so large that we can assume it 
to be unaffected by any of the reactions we consider, thus 
we fix pion yield in terms of fugacity and temperature 



values. We do not consider all resonances which decay 
into resonances as does the SIIARE2 SHM program \Q\, 
thus we will correct the final yields by an estimate of this 
effect comparing our initial resonance yield to SHARE2 
results. 

For the 'fast' baryon resonances considered here we 
keep the sum yields constant: 
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The baryon annihilation, strangeness exchange such as 
N+K <-^- A-|-7r reactions, and population exchanges with 
higher resonances are assumed not to have a material im- 
pact within the time scale during which the temperature 
drops from chemical to kinetic freeze-out condition. 

The experimentally observable hyperon yield appear- 
ing in our final result is 



Atot = S(1385) + A + S°(1193) + Y* 



(2) 



due to experimentally inseparable S°(1193) — ^ 7 + A 
decay and the decay of further hyperon resonances Y*. 
Similarly, when we refer to A^tot we include baryon reso- 
nances in the count. 

In the following we will be referring explicitly to the 
A yield governed by cta = I/Fa = 1.67 fm. All equa- 
tions apply equally to E(1385) yield (partial decay width 
Fs^A — 35 MeV) and we will compare our results with 
experiment for this case. We note that even though the 
S(1385) decay width is much smaller than Fa, the num- 
ber of reaction channels and particle densities available 
lead to a significant effect for S(1385), comparable to our 
finding for A. 

The evolution in time of the A (or £(1385)) resonance 
yield is described by the process of resonance formation 
in scattering, less natural decay: 
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where A^a is multiplicity of A resonances, 
dWN-K^A/dVdt and dWA-^N-n/dVdt are invariant 
rates (per unit volume and time) for A production and 
decay respectively. Allowing for Fermi-blocking and 



Bose enhancement in the final state, the two in-niatter 
rates are: 
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where g^, i = tt, A^, A is particles degeneracy. The distri- 
bution functions for tt, iV, A are 



h = 



1 



T-ie«-P-/'r-l' 
1 



T7ie«-P^/T' + i 



, J = A^,A. 



(6) 
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Here T^ is particles fugacity, and u ■ pi — Ei, for 
u^ = (1,0) in the rest frame of the heat bath where 
d^p5o(pf — "T-f) — > d^pi/Ei for each particle. Hence, 
Eq.(|3]) and Eq.® are Lorentz invariant, as presented 
these rates can be evaluated in any convenient frame of 
reference. Normally, this is the frame co-moving with the 
thermal volume element. 

Eqs.®, ([7]) satisfy the usual relations for the Fermi 
and Bose distributions: 

1-f, - Tj'e^-P^f,, j = A,7V (8) 

i + U = r-'e'-^'U- (9) 

Using these equations, we obtain for decay rate Eq.(|3]) 
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where in the last exponent we replaced pm + P-k by pa 
given the energy-momentum conservation 4-delta func- 
tion. We can perform a similar simplification in Eq.(l5]). 
Then, observing that due to the time reversal symmetry, 

|(pi|Mb2P3)|' = |(p2P3|A/bi)|' 



we find a detailed balance relation between the produc- 
tion and decay rate: 
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The master equation, Eq.Q, can now be cast into the 
form: 
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This is a rather intuitive and simple result, yet only re- 
cently the 1^ 2 population master equations have been 
considered ill]. Equation P^ implies for dNA/dt = 
the chemical equilibrium condition: 



^TT ^ N ~ ^ A- 



(13) 



This equation is solved by the global chemical equilib- 
rium T'^^ — T^^ — T^ = 1. However, there are also 
other, transient, equilibrium states possible, given a pre- 
scribed value of e.g. the background pion abundance, 
T'^ y^ 1. When the initial state is formed away from 
transient equilibrium condition, we recognize that for 
Ta < TttTjv the A production is dominant, and con- 
versely, for Ta > TttTat the A decay dominates. 

We now introduce into the population master equation 
(|12p the effective lifespan, ta aiming to find an equation 
similar to classic radioactive decay population equation. 
We define the in medium A-lifespan to be: 



TA 



Ta dNA/dTA 
V dWA-,N^/dVdt' 



(14) 



We recognize that in the Boltzmann limit this corre- 
sponds to the ratio of equilibrium yield to the rate per 
unit time at which the equilibrium is approached. We 
obtain for Eq.ini): 
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In case that the ambient temperature does not vary 
with time, and thus only populations evolve due to 
change in fugacitics, we have dN /dt ~ dN/dT dT /dt and 
the following dynamical equation for the fugacity arises: 
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This is 'classical' population equation form where the 
fugacity plays the role of the classical densities. When 
the dynamical values of Ti(i) are used in the quantum 
Bose/Fermi distributions, the effects of blocking, and 
stimulated emission are explicit. 

If we instead were to introduce the lifespan by ta = 
{N a/V) / [dWA^N-n / dV dt) , this imphes for all particles 
(Bose, Fermi, Boltzmann) the classical population equa- 
tion, e.g. dNA/dt — (T^Tat/Ta — 1)A^a/ta, and the 
quantum effects are now hidden in the definition of r. 



Both definitions coincide for the case of a dilute system, 
and differ most for dense systems. In the limit of very 
dilute, vacuum system, the relaxation time is the same 
as the lifespan of the particles. The computed yields of 
particles as function of time are not dependent on the 
finesse of the relaxation time definition. 

We now set up for semi-analytical solution of master 
equation ^TE\i . For multiplicities A and N considering 
the small yield and m ^ T we will use the Boltzmann 
distribution: 
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(17) 
(18) 



where xa,n — f^AjN/T, K2{x) is Bessel function. Con- 
sidering that fugacities, temperature and volume vary in 
time, we rewrite the left hand side of Eq. (fT5|) : 
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We changed from t to r to make explicit the fact that 
we work in fluid-element co-moving frame and thus do 
not consider the effect of flow on the volume time depen- 
dence. 

Combining Eq. p3|) with Eq. p^ we obtain 
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The last term is negligible, T5 ^ ta, tt since pions dom- 
inate and we have near conservation of entropy which for 
massless particles would in fact imply VT^ = Const. 

Since entropy must be (slightly) increasing, while T is 
decreasing with time, rg > 0. Similarly, tt > 0, since the 
temperature decreases with time, and x'^K2(x), x — rajT 
increases with T : 



x^K2{x) w Vo^x^^"^ exp{-x); 



Therefore: 
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We now evaluate the magnitude of tt invoking a model 
of matter expansion of the type used e.g. in [13|], where 
the longitudinal and transverse expansion is considered 
to be independent. In the proper rest frame of the out- 
flowing matter. 



fo.T^^=nRl{T)T'^ 
dy dy dy 



Const.. 



(25) 



We will use dz/dy ~ r, where r is the proper time in 
the local volume element, this is exact for a 1-d ideal 
hydro flow. The growth of the transverse dimension can 
be generically described by: 



R±{t)=Ro+ [ 

J Tn 



v{T')dT'. 



From Ea. (|25p by elementary evaluation we obtain: 

t _ 1 (2{vt/Ri_) + 1 
T ~ ^3 



(26) 



(27) 



Equation (j27p evaluated near hadronization condition 
is yielding the magnitude of tt, see Eq. ([M|) . If the maxi- 
mum expansion velocity is practically instantly achieved, 
vt/R± ~ 1. This leads to maximum value of T/T ~ 
— 1/r. However if a more realistic proflles are assumed, 
T/T is diminished in magnitude as much as 30%. We 
thus conclude that 

0.5 TTiA 1 0.7 toa 

Th T TT Tfi T 

which for hadronization time t^ < 10 fm can compete 
with the width of the A-resonance, I/^a — 120 MeV. As 
this shows, the details of the expansion model are not 
critical for the results we obtain. In actual calculations 
we employ w(r) described in [l3|, where we assume that 
the expansion is already at maximum velocity at the time 
of chemical freeze-out. The resulting dependence T(t) 
after chemical freeze-out is shown in figure [TJ We note 
that the time between chemical and thermal freeze-out 
At is not longer than about 2.5fm/c, and can be as short 
as Ifm/c. However, even such a short scattering period 
is enough to alter the visible yields of strong resonances, 
in fact most pronounced effect we find in the latter case, 
since the longer time allows a greater degree of chemical 
equilibration. 

We now can solve Eq. ipO)) . Employing Eq.([T|) we have: 
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(29) 



(30) 



where N^ and N^ are densities of A and N resonances 
with Ta = Tjv = 1. The solution of Eq. ipS)) is elemen- 
tary: 



Ta(t)= T° 



•'^'- dT I e ■'^h (31) 



qe-'^h 



where t/,, is initial expansion time at hadronization, and 
Th < T < Tniax, Upper time limit chosen to yield Tmax = 
120 MeV, i.e. r^ax ^ 8 fm. 
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FIG. 1: Temperature T as function of (5r, the proper 
time interval between chemical and thermal freeze-out or 
chemical freeze-out temperature (from top to bottom) T = 
180, 160, 140 MeV and thermal freeze-out T > 120 MeV. 



FIG. 2: The ratio A/A° (solid lines) and iV/iV° (dashed hues) 
as functions of temperature T for select given pairs of values 
T,'yq, see text and figure box for details. 



In order to evaluate the final A multiplicity we 
need also to know initial particles densities right after 
hadronization which we consider for RHIC head-on Au- 
Au collisions at ^snn = 200 GeV. We introduce the 
initial hadron yields inspired by a picture of a rapid 
hadronization of QGP with all hadrons produced with 
yields governed by entropy and strangeness content of 
QGP by quark recombination. In this model the yields 
of mesons and baryons are controlled by the constituent 
quark fugacity 7^: 









(32) 



Thus for 7q > 1 we have the condition Ta < TttTat, and 
the yield of A will increase in the time evolution. 

For each entropy content of the QGP fireball, the corre- 
sponding fixed background value of 7^ can be found once 
hadronization temperature is known [12| . For T = 140 
MeV pions form a nearly fully degenerate Bose gas 
with 7q ~ 1.6. In the following discussion, aside of 
this initial condition, we also consider the value pairs 
T = 150 MeV, 7, = 1.42, T = 160 MeV, 7, = 1.27, 
T = 170 MeV, 7q = 1.12 and T = 180 MeV with 7, = 1. 
Note further that for m ^ T the density A is relatively 
low, thus there is no significant dependence of I/ta on T 
and Ta; in essence ta = /i/Fa takes the free space value 
TA ^ ;i/120MeV. 

In figure [2] we present results for ratios A/Aq (solid 
lines) and N/Nq (dashed lines) as functions of tempera- 
ture T, beginning from the presumed initial hadroniza- 
tion temperature T through Tmax = 120 MeV. Aq and 
A'o are the initial yields obtained at each hadronization 



temperature. For T < 180 MeV, initially Ta < TatT^t, 
thus based on our prior discussion, we expect that the 
master equation leads to an initial increase in the yield 
of resonances. However, as temperature drops, due to 
the dynamics of the expansion the increasing yield of A 
turns over, and a final nett increase of resonance yield is 
observed for T < 160 MeV. We note that for T > 180 
MeV there is a continuous depletion of resonance yield. 
The nucleon yields move in opposite direction to the A- 
resonance. 

This behavior can be understood in qualitative man- 
ner as follows: The total number A -I- A^ is conserved 
therefore A multiplicity increases and A^ multiplicity de- 
creases until they reach transient chemical equilibrium 
{dN/\/dT = 0), corresponding to the maximum point 
seen for A in figure [D There is also infiuence of expan- 
sion: even if for some temperature the transient equilib- 
rium condition ()13p is reached, the system cannot stay in 
this equilibrium, Ta and Tjy are increasing to conserve 
total number of particles. Ta increases faster because of 
larger A's mass. After Ta becomes larger than TttTat 
A decay begins to dominate and their multiplicity is de- 
creasing. The special case at hadronization temperature 
T — 180 MeV where, T^ = 1 and equilibrium condition 
is satisfied initially. As expansion sets in, A is decreas- 
ing because Ta > Tjv (recall that here T^r = !)■ In the 
SHM evaluation of yields one assumes that all ratios seen 
in figure [2] are unity. 

The initial hadronization yields which we used as ref- 
erence in figure [2] are not accessible to measurement. 
Therefore, we consider in figure [3] the fractional yield 
A/Atot (top frame), again as a function of temperature 
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FIG. 3: Relative resonance yield, for (top) A/A'tot and (bot- 
tom) E(1385)/Atot as a functions of freeze-out temperature, 
for hadronization temperatures To = 140, 160, 180 MeV, see 
box and text for details. The dotted brown line gives the 
expected SHM chemical equilibrium result. 



T. The results for hadronization temperatures Tq = 140 
(solid blue line), Tq — 160 (dash-dot green line) and 
To — 180 MeV (dashed brown line) are shown. TVtot is 
fixed by hadronization condition and is not a function of 
time, as discussed. Thus the observable final rapidity nu- 
cleon yield corresponds to the initial value at hadroniza- 
tion. Note that up to strange and multi strange baryon 
contribution, A'tot is the total baryon (rapidity) yield. 

Since in this study we have considered a subset of 
all relevant baryon resonances our chemical equilibrium 
reference yield (line for Ta = Tn) is not the same 



as the corresponding reference line for the full statisti- 
cal hadronization model (SHM) evaluation, obtained us- 
ing SHARE2, and presented as 2A+/p (upper frame) 
and 3E+/A0 (lower frame). The SHARE2-SHM value 
A++/p ~ 0.2 at T ~ 160 MeV is consistent with the 
STAR d-Au results [l3|. Also, comparing our with the 
SHARE2 result we note that SHARE2 yield is larger at 
chemical freeze-out. The magnitude of the difference in 
the yields at time of chemical freeze-out provides a mea- 
sure of the magnitude of the corrections we can expect to 
arise in the full treatment at thermal freeze-out and/or 
systematic error for these yields. 

The nature of these effects is different for the two 
yield cases considered: the presence of heavier resonances 
which cascade by way of A leads to an increase of the 
thermal freeze-out yield. The correction is thus nearly 
as much as we see the SHARE2 yield higher at chemi- 
cal freeze-out. For E(1385) the difference with SHARE2 
arises from a difference of contributions of partial decays 
producing Atot, thus the correction is multiplicative fac- 
tor which does not change, but is uncertain in magnitude 
due to lack of knowledge about the branching ratios. 

We believe that the A and S(1385) yields are un- 
derestimated by about 15% - 35%. (bigger effect for 
hadronization at higher T). This implies that depend- 
ing on hadronization temperature a relative yield range 
0.16 < A++/p = 0.5A/7Vtot < 0.26 arises, and similarly 
(see lower frame in figureED 0.35 < S](1385)/Atot < 0.43 
with the higher relative yield corresponding to the lower 
hadronization temperature. One of the key results of 
this work is the narrow range for E(1385)/Atot, and the 
fact that the initial chemical non-equilibrium effect leads 
to a reversal of the SHM model situation: the relative 
yields of massive resonances decreases with decreasing 
hadronization temperature. . 

In order to compare with the experimental results 
we note that the data presented [l|, 0| are for charged 
S(1385), particle and antiparticle channels, (S^(1385)-|- 

E±(1385)/(Atot +M^t) ^ 0.29. This result needs to be 
multiplied with 3/2 to be comparable to results presented 
here which include I]°(1385). Multiplying the value for 
hadronization at T = 140 MeV with thermal freeze-out 
a.t T — 120 MeV, and allowing for contribution by heav- 
ier resonances as indicated by SHARE2 our result is in 
perfect agreement with [1, 2] However, given the narrow 
range of results we find, it seems that the high yield of 
11(1385), seen the error 0(20%) is nearly compatible with 
the entire range of chemical freeze-out temperatures here 
considered - the low T chemical freeze-out is favored by 
1.5 s.d. over high T. 

The reader should take note that the 'thermal' model 
result presented in Ref. !l| corresponds to initial high 
temperature freeze-out in chemical equilibrium which is 
unobservable, since the high T hadronization resonance 
decay products have no chance to escape into free space. 
Thus this comparison of this model with experiment is 
flawed. The evolved yield is shown as (red) dashed line in 
figureO and is found 25% below the value measured. The 



reason this happens is that the high T chemical freeze- 
out happens near chemical equilibrium and the yields 
follow closely the chemical equilibrium yield described by 
temperature, thus it is the thermal freeze out temperature 
which in this case controls the final observable resonance 
yield. 

In summary, we have presented master equation gov- 
erning the evolution in time of the A, I](1385) baryon 
resonance yield after QGP hadronization, allowing for 
resonance decay and production process. We have shown 
considering the properties of the master equation that if 
the yield of hadrons is initially above chemical equilib- 
rium, the resonance population increases beyond the ini- 
tial yield. Conversely, we find that in a physical system 
in which the particle multiplicities of hadrons arise be- 
low chemical equilibrium yields, a circumstance expected 
below threshold to QGP formation, the final yield of res- 
onances is suppressed by the dominance of the resonance 
decay process over back reaction resonance production. 

In a quantitative model we evolved the yields after 
QGP hadronization allowing for initial chemical non- 
equilibrium particle abundances, and volume expan- 
sion assuring entropy conservation. We found, see fig- 
ure [21 that the thermal freeze-out fractional resonance 
yield differs significantly from the chemical-freeze out 
SHM expectation, with the scenario involving high-T 
hadronization resonance yield being depleted, and low-T 
hadronization yield scenario further enhanced in relation 
to the total yield. 

The resonance enhancement effect we presented can 
only occur when the initial state is out of chemical 
equilibrium, and the decay/formation processes are fast 
enough to compete with the hadron volume evolution. 
One would thus think that 'narrow', i.e. quasi-stable 
resonances are not subject to the effects considered here. 
However, a special consideration must be given to nar- 



row resonance which are strongly coupled to more mas- 
sive resonances which can decay fast into other chan- 
nels. An example is A(1520). In this situation the lower 
mass state, here A(1520), is the ground state which is de- 
pleted by coupling to the yet more massive state. There 
are several states of relevance to consider, hence A(1520) 
'quenching' 7] must be addressed in a more rigorous nu- 
merical approach. Wc will to return to discuss this case 
in a separate publication [l5|. 

Aside of several specific prediction we made here, there 
are three important general consequences of our study: 
a) the fractional yield of resonances A* /A can be con- 
siderably higher than expected naively in SHM model of 
QGP hadronization, b) since there is nearly a factor two 
difference in the final thermal freeze-out ratio in A/A'tot , 
while the SHM yields a more T independent result, one 
can imagine the use of A/A^tot as a tool to distinguish 
the different hadronization conditions e.g. chemical non- 
equilibrium vs chemical equilibrium a point noted in sim- 
ilar context before [ig| ; and c) we have shown that the rel- 
atively high yield of charged S* (1385) reported by STAR 
is well explained by our considerations with hadroniza- 
tion at T = 140 MeV being favored. 
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